Fermionic Hopf solitons and Berry's phase in topological surface superconductors 
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A fascinating idea in many body physics is that quantum statistics may be an emergent property. 
This was first noted in the Skyrme model of nuclear matter, where a theory formulated entirely 
in terms of a bosonic order parameter field contains fermionic excitations. These excitations are 
smooth field textures, and believed to describe neutrons and protons. We argue that a similar 
phenomenon occurs in topological insulators when superconductivity gaps out their surface states. 
Here, a smooth texture is naturally described by a three component real vector. Two components 
describe superconductivity, while the third captures the band topology. Such a vector field can as- 
sume a 'knotted' configuration in three dimensional space - the Hopf texture - that cannot smoothly 
be unwound. Here we show that the Hopf texture is a fermion. To describe the resulting state, the 
regular Landau-Ginzburg theory of superconductivity must be augmented by a topological Berry 
phase term. When the Hopf texture is the cheapest fermionic excitation, interesting consequences 
for tunneling experiments are predicted. 
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There has been much recent excitement relating to 
topological insulators (TIs), a new phase of matter with 
protected surface states P, Q ■ Particularly rich phenom- 
ena are predicted to arise when this phase is combined 
with conventional orders such as magnetism 0, 3] , crys- 
talline order [HI and superconductivity. The last is par- 
ticularly interesting. Superconductivity induced on the 
surface of a TI was predicted to have vortices harboring 
Majorana zero modes 0, similar to the bound states in 
vortices of a p + ip superconductor 0-0] . These are of 
interest to quantum information processing, since they 
are intrinsically robust against errors. Receritly, super- 
conductivity was discovered in a doped TI 'iq, which 
could be used to induce surface superconductivity. Be- 
low we discuss a new theoretical approach to studying 
this remarkable superconducting phase, which provides 
different insights and directions for experiments. 

We focus on smooth configurations where the energy 
gap never vanishes. In this case, the low energy descrip- 
tion of the system is entirely in terms of bosonic 'order 
parameter' coordinates, much as the Landau Ginzburg 
order parameter theory describes superconductors at en- 
ergies below the gap. Can fermions ever emerge is such 
a theory? While it is easy to imagine obtaining bosons 
from a fermionic theory, the reverse appears impossible at 
first sight. However, it has been shown in principle that 
bosonic theories that contain additional Berry's phase (or 
Wess-Zuniino-Witten) terms, can accomplish this trans- 
mutation of statistics. We show that this indeed occurs 
in the superconductor- TI (Sc-TI)system; the order pa- 
rameter theory contains a Berry phase term which im- 
plies that a particular configuration of fields - the Hopf 
soliton (or Hopfion) - carries fermionic statistics. While 
such statistics transmutation is common in one dimen- 
sion eg. bosons with hardcore interactions 
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matter context, an physically realizable example exists 
in two dimensions: solitons of quantum Hall ferromag- 
nets (skyrmions) are fermionic and charged, and have 
been observed [l2|-[l4[. However, the superconductor-TI 
system is, to our knowledge, the first explicit condensed 
matter realization of this phenomenon in three dimen- 
sions. 



MODEL AND HOPF TEXTURE 

The essential properties of a topological insulator are 
captured by a simplified low energy theory with a three 
dimensional Dirac dispersion (a microscopic realization 
is described later): Hd = ip'^lvpct ■ p -I- m/?o]V'! where 
(ai, a2, Q!3, /3o) are 4x4 anti-commuting matrices which 
involve both spin and sublattice degrees of freedom. The 
matrices ai can be taken to be symmetric, while /?o and 
/?5 — Q;iQ;2a3/?o are antisymmetric. The dispersion then 
is e{p) 



±yjvpp^^~\-rn? . An insulator is obtained for 
TO 7^ 0. Changing the sign of to results in going from 
a trivial to a topological insulator. Which sign of 'm' is 
topological is set by the band structure far away from the 
node - here we assume to < is topological. Since the 
vacuum can be taken to be a trivial insulator, this mass 
term changes sign at the topological insulator surface. 
Consider now adding (onsite) superconducting pairing, 
which may be proximity induced by an s-wave supercon- 
ductor. Then iJpair — A7/'t/35'0t -I- h.c, and we can write 
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rare phenomenon in higher dimensions. In the condensed 



(1) 

the spectrum now is e(p) = ± v^p'^ + where 
A'P = m? + |Ap. It is convenient to define a three vec- 
tor n = (T^gAjItoA, to), such that \n\ — M. A singular 
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FIG. 1: The Hopf map: fn : r — > n is shown by displaying 
contours of equal h. Points at infinity are all mapped to the 
same point on the sphere fnioo) = z. In red is fn^ln = z], 
in green fjj^[—z], in blue f^^[x], and in yellow f^^ly]- Note 
the unit linking of any pair of curves, which can be used to 
define the Hopf texture. 

configuration is one where all three components of this 
vector go to zero. This corresponds to the core of a vortex 
(where the components of the pairing A vanish), inter- 
secting the surface of a topological insulator (where the 
third component is zero). It has been pointed out that 
an odd strength vortex will give rise to an unpaired Ma- 
jorana zero mode in this configuration [6^ , which can also 
be viewed as a hedgehog with odd integer topologi- 
cal charge. Note, at the core of these singular configura- 
tions the gap closes, allowing for the possibility of local- 
ized bound states near zero energy. In this work we will 
only consider smooth textures of the n field, where the 
single particle gap is nonzero everywhere. An effective 
theory of slow fluctuations of the 'order parameter' field 
n{r, t) (occuring over spatial (time) scales much larger 
than ^ = hvp/M (r = £,/vf)), can be obtained by in- 
tegrating out the gapped fermions. An analogous proce- 
dure is well know in the context of the BCS theory of su- 
perconductivity, where it leads to the Landau- Ginzburg 
action. Here we will find that an extra topological term 
arises, that transmutes statistics and leads to fermionic 
solitons. 

Consider a smooth configuration of n(r) , which can be 
normalized to give a unit vector n(r) at each point. This 
defines a mapping from each point of three dimensional 
space, to a unit three vector, which describes the surface 
of a sphere . We require that the mapping approaches 
a constant at infinity: n{\r\ oo) = const (e.g the vac- 
uum) . Can all such mappings be smoothly distorted into 
one another? A surprising result due to Hopf [3| 1931, 
is that there are topologically distinct mappings, which 
can be labeled by distinct integers h (the Hopf index) . No 



smooth deformation can connect configurations with dif- 
ferent Hopf indices. Mathematically, Hopf showed that 
the homotopy group: H3[S'^] ~ Z. A straightforward 
way to establish the index is to consider the set of points 
in space that map to a particular orientation of h. In 
general this is a curve. If we consider two such orienta- 
tions, we get a pair of curves. The linking number of the 
curves is the Hopf index. A configuration with unit Hopf 
index can be constructed by picking a reference vector 
n — z, say, and rotating it by the following set of rota- 
tions. Any rotation is parameterized by an angle and a 
direction of rotation. If we take the angle to vary as we 
move in the radial direction, from 6* = at the origin, 
to 9 = 211 at radial infinity, and take the axis of rotation 
to be the radial direction f, this gives a Hopf texture. 
Using an SU(2) matrix U{'r) = e^T-'" '^ to represent this 
rotation, the vector field n(r) — U {r)azU'^ if) is the unit 
Hopf texture. This is readily verified by studying which 
spatial points map to fi = ±z. While the former includes 
all points at infinity as well as the z axis, the latter is a 
circle in the xy plane. Clearly these two curves have unit 
linking number unity (see Figure [TJ What is the physical 
interpretation of this Hopf texture in the context of TIs? 
A torus of TI (Figure [2|) has superconductivity induced 
on its surface. There is vacuum far away and through 
the hole of the torus, which counts as a trivial insulator, 
n = z. The center of the strong topological insulator cor- 
responds to n = —z. On the topological insulator surface 
= 0, and the superconducting phase varies such that 
there is a unit vortex trapped in each cycle of the torus. 
We now argue that such a texture is a fermion. 

HOPF SOLITONS ARE FERMIONS 

The ground state of the mean field Hamiltonian ([T]) for 
a general texture, has rather low symmetry, and cannot 
be labeled by spin or U(l) charge quantum numbers. The 
only quantum number that can be assigned is the parity 
of the total number of fermions (—1)^. Superconduct- 
ing pairing only changes the number of fermions by an 
even number, hence one can assign this Z2 fermion par- 
ity quantum number to any eigenstate. We now argue 
that the fermion parity of a smooth texture is simply the 
parity of its Hopf index. 

First, we argue that the ground state with a topologi- 
cally trivial texture has an even number of fermions. As 
a representative, consider a configuration where the su- 
perconductor pairing amplitude is real. This is a time 
reversal invariant Hamiltonian. If the ground state had 
an odd number of fermions, it must be doubly degenerate 
at least, by Kramers theorem. However, the ground state 
of any smooth texture is fully gapped and hence unique. 
Thus, this configuration must have an even fermion par- 
ity. Now, any other texture in the same topological 
class can be reached by a continuous deformation, dur- 



FIG. 2; This figure shows a torus of strong topological insu- 
lator (green) in vacuum, whose surface is superconducting . 
There is unit superconductor phase winding about each cycle 
of the torus. We plot the equal phase contours on the sur- 
face whose pairing phase is O(blue) and 7r(yellow). The unit 
linking of these curves indicates this is the Hopf mapping. 

ing which the gap stays open. The fermion parity stays 
fixed during this process. Note, this argument cannot be 
applied to configurations with nonzero Hopf index, since 
these necessarily break time reversal symmetry. For ex- 
ample, the configuration shown in Figure 1, contains a 
pair of vortices. 

To find the fermion parity of the nontrivial Hopf config- 
urations, we consider evolving the Hamiltonian between 
the trivial and h = 1 configuration. In this process we 
must have n = at some point, which will allow for the 
gap to close, and a transfer of fermion parity to poten- 
tially occur. Indeed, as shown below in separate calcu- 
lations, a change in fermion parity is induced when the 
Hopf index changes by one. 

Numerical Calculation: We study numerically the mi- 
croscopic topological insulator model defined in refer- 
ence JTJ, with a pair of orbitals (t^ = ±1) on each 
site of a cubic lattice. The tight binding Hamiltonian 
^ ^ Sfcl'/'fc^feV'fc + -f^pair(^)], is written in momen- 
tum space using a four component fermion operator 
Tpk with two orbital and two spin components. Then, 
Hk = -2tYfa=i^a sin ka - mlSolX + J2l^iicoska ~ 1)], 
where (a, /3o) = [txCTx, r^ay, -Ty.T^). For m > a 
strong topological insulator is obtained when A £ (0,2). 
In addition we introduce onsite singlet pairing ^^pair(^) — 
^[V'il'^craV'lfc + h.c. Note, when A « 0, fc « (0, 0, 0), 
Eqn. [1] is recovered as the low energy theory. 

The energy spectrum is studied as we interpolate be- 
tween a topologically trivial texture {h = 0) and the Hopf 
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FIG. 3: The spectral flow of the lowest 20 eigenvalues when 
the pairing on the surface of the topological insulator in Fig. 
[2] is linearly interpolated between two limits: A(a:;, A) = (1 — 
A)Ao(a;) -I- AAi(a;). Ao(x) is constant over the whole surface. 
TOP: Ai(i) has a unit phase winding (vortex) in both the Ri 
and i?2 cycles of the torus, i.e. the Hopf texture. BOTTOM: 
Ai(a;) has a unit phase winding in only the R\ cycle. It is 
clear that there is a single level crossing in the TOP case, 
meaning the ground state fermion parity is changed in the 
process. The initial state has even fermion parity, so the final 
state, the Hopf texture, must have odd fermion parity. The 
calculation is on the cubic lattice model defined in the text, 
with R\ — 13 and f?2 = 5 lattice units, and only points within 
the torus are retained. Parameters used: t = M = X = 1, and 
pairing |A| = 1. 



texture {h = 1). We choose to define a torus shaped 
strong topological insulator with trivial insulator (vac- 
uum) on the outside, as in figure[2j The surface is gapped 
by superconducting pairing A, which in the trivial tex- 
ture is taken to be real Aq. In the Hopf texture, the su- 
perconducting pairing Ai has a phase that winds around 
the surface, with a unit winding about both cycles of the 
torus. This can be interpreted as 'vortices' inside the 
holes of the torus. Note, the vortex cores are deep in- 
side the insulators, so there is a finite gap in the Hopf 
texture. We interpolate between these two fully gapped 
phases be defining A(A) — AAi + (1 — A) Aq and changing 
A = — > 1. On the way, the gap must close since the 
two textures differ in topology. We study the evolution of 
eigenvalues as shown in figure[3l We find that exactly one 
pair of ±E eigenvalues are pumped through zero energy. 
As argued below, this signals a change in the fermion 
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parity of the ground state on the two sides. Since the 
trivial texture has even fermion parity from time rever- 
sal symmetry, the Hopf texture must carry odd fermion 
number. No such odd level crossings occur for topolog- 
ically trivial textures (eg. phase winding through only 
one cycle of the torus) . 

To see why the crossing of a ±E conjugate pair of levels 
corresponds to a change in fermion number, consider a 



single site model H = Eo{c'c 



This has a pair of 



single particle levels at ±Eo, which will cross if we tune 
Eq from say positive to negative values. However, writing 
this Hamiltonian in terms of the number operators H = 
Eo{2n — 1), shows that the ground state fermion number 
changes from n = to n = 1 in this process. Thus the 
ground state fermion parity is changed whenever a pair 
of conjugate levels cross zero energy. 

The Pfajjian: Previously, the ground state fermion 
parity was found by interpolating between two topologi- 
cal sectors. Can one directly calculate the fermion parity 
for a given Hamiltonian's ground state? We show this is 
achieved by calculating the PfafHan of the Hamiltonian 
in the Majorana basis. The Pfaffian of an antisymmet- 
ric matrix is the square root of the determinant - but 
with a fixed sign. It is convenient to recast the Hamil- 
tonian in terms of Majorana or real fermions defined via 
V'a = (xia + *X2a)/2. Siuce a pair of Majorana fermions 
anticommute {xi, Xj} — ^iji the Hamiltonian written in 
these variables will take the form: 



H 



XiXj 



(2) 



where is an even dimensional antisymmetric matrix, 
with real entries and the Majorana fields appear as a vec- 
tor X — {■ ■ ■ 1 Xia, X2a, ■ • • ); whcrc a refers to site, orbital 
and spin indices. The iE symmetry of the spectrum is 
an obvious consequence of h being an antisymmetric ma- 
trix. The ground state fermion parity in this basis is 
determined via: 



(-1)^ = sign [Pfaffian(h)] 



(3) 



We numerically calculated the Pfaffian of a Hamilto- 
nian with a single Hopf texture for small systems and 
confirmed it has a negative sign. In contrast, the triv- 
ial Hopf texture Hamiltonian has positive Pfaffian in the 
same basis. 

Finally, we mention that it is possible to confirm the 
numerical results analytically, by solving for the low en- 
ergy modes in the vicinity of the vortex core, where the 
insulating mass term is set to be near zero. The linking 
of vortices in the Hopf texture plays a crucial role in de- 
riving this result. In the Appendix, we discuss how this 
result is connected to the three dimensional non-Abelions 
in the hedgehog cores of Ref [ll. 

A Two Dimensional Analog: We briefly mention a two 
dimensional analog of the physics described earlier. Note, 



eqn. [T] with the third component of momentum absent 
Pz = 0, describes a quantum spin Hall insulator (triv- 
ial insulator) when m > (to < 0), in the presence of 
singlet pairing A. Again, as before a three vector charac- 
terizes a fully gapped state, and the nontrivial textures 
are called skyrmions (H2(S'2) = Z). A unit skyrmion 
can be realized with a disc of quantum spin Hall insu- 
lator with superconductivity on the edge, whose phase 
winds by 2tt on circling the disc. Again, one can show 
that the skyrmion charge Q determines the fermion par- 
ity (—1)*^ = (— 1)'^. An important distinction from the 
three dimensional case is that the low energy theory here 
has a conserved charge. If instead of superconductiv- 
ity, one gapped out the edge states with a time reversal 
symmetry breaking perturbation which had a winding, 
then this charge is the electrical charge. It is readily 
shown that the charge is locked to the skyrmion charge 
Q. Hence odd strength skyrmions are fermions [s^. This 
is closely analogous to the Quantum Hall ferromagnet, 
where charged skyrmions also occur Returning 
to the case with pairing, since that occurs on a one di- 
mensional edge, it is difficult to draw a clear-cut separa- 
tion between fermions and collective bosonic coordinates, 
in contrast to the higher dimensional version. Hence we 
focus on the 3DTIs. 



EFFECTIVE THEORY AND TOPOLOGICAL 
TERM 

The gap to the iJj fermions never vanishes since |n| > 0, 
so one can integrate out the fermions to obtain a low 
energy theory written solely in terms of the bosonic or- 
der parameter n. How can this field theory describe a 
fermionic texture? As described below, this is accom- 
plished by a topological Berry phase term which appears 
in the effective action for the n field. 

In computing the topological term, it is sufficient to 
consider a gap whose magnitude is constant r?(r, t) — 
Mn{r,t). Integrating out the fermion fields with action 
Sf = J d* x[ip^dtip — Hf], (the integral is over space and 
(Euclidean) time) , one obtains the effective action for the 
bosonic fields: 



-SBiii{r.t) _ 



(4) 



This computation may be performed using a gradient 
expansion, i.e. assuming slow variation of the n field 
over a scale set by the gap. Two terms are obtained: 
Sb = So + Stop- The first [3l| is a regular term that pe- 



nalizes spatial variation: Sq = /(9^n)^. The second 
is a topological term which assigns a different amplitude 
to topologically distinct spacetime configurations of h. 
Assuming ri(oo) = const, these configurations are char- 
acterized by a Z2 distinction (H4(S'2) = Z2)[18]. That 
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is, there are two classes of maps - the trivial map, which 
essentially corresponds to the uniform configuration, and 
a non-trivial class of maps, which can all be smoothly 
related to a single representative configuration ni{r,t). 
If the function r[{n(r, t)}] = 0, 1 measures the topolog- 
ical class of a spacetime configuration, then the general 
form of the topological term is 5'top = i9T[{n{r, t)}]. The 
topological angle 9 can be argued to take on only two 
possible values 0, tt, since composing a pair of nontrivial 
maps, leads to the trivial map. Via an explicit calcula- 
tion, outlined below, we find 6 = n. Let us first examine 
the consequences of such a term. The nontrivial tex- 
ture ni(r, can be described as a Hopfion-antiHopfion 
pair being created at time ti, the Hopfion being rotated 
slowly by 2tt, and then being combined back with the 
anitHopfion at a later time t2 [19]. The topological term 
assigns a phase of to this configuration. This is equiv- 
alent to saying the Hopfion is a fermion, since it changes 
sign on 27r rotation. 

Calculating the topological term requires connecting 
the pair of topologically distinct configurations. To do 
this in a smooth way keeping the gap open at all times 
requires enlarging the order parameter space for this pur- 
pose. If m(r, t, A) is an element of this enlarged state that 
smoothly interpolates between the trivial configuration 
m(r, t, 0) — const, and the nontrivial one m(r, i, 1) — 
rii(r, t), as we vary A, then one can analytically calcu- 
late the change in the topological term dStop/dX, and 
integrate it to get the required result [IB-Iil. The key 
technical point is finding a suitable enlargement of our 
order parameter space 82- Remembering that this can 
be considered as S2 = SU{2)/SO{2), we can make a 
natural generalization A/3 = SU{3)/SO{3). The latter 
has all the desirable properties of an expanded space, eg. 
there are no nontrivial spacetime configurations, so ev- 
erything can be smoothly connected (114 (M3) — 0). This 
extension allows us to calculate the topological term, (as 
explained in detail in the methods section), which yields 



PHYSICAL CONSEQUENCES 

We now discuss two kinds of physical consequences 
arising from fermionic Hopfions. The first relies on the 
dynamical nature of the superconducting order parame- 
ter, while the second utilizes the Josephson effect to iso- 
late an anomalous response. The first class conceptually 
parallels experiments used to identify skyrmions in quan- 
tum hall ferromagnets. There, when skyrmions are the 
cheapest charge excitations, they are detected on adding 
electrons to the system 

Consider surface superconductivity on a mesoscopic 
torus shaped topological insulator as in Figure [2] The 
Hopf texture corresponds to unit phase winding in 
each cycle of the torus. The energy cost. Eh = 




$1 — normal 



2-K 



/47r 



27r— = TT anomalous 




odd fermion even fermion 
parity parity 



FIG. 4: Anomalous Josephson Response connected to 
fermionic Hopfions: An annular cylinder of TI with pairing 
induced on the top and bottom surfaces via proximity to a 
superconductor (present in the red regions) is shown. Tun- 
ing the flux to $1 ~ V'o/S, induces a vortex in the Ci cycle. 
Now, the Josephson effect on tuning $2 will be anomalous, 
with a part that is not periodic in the flux ipo. This is di- 
rectly related to the change in ground state fermion parity at 
$2 = 27r((po/27r), where the Hopf texture is realized. Adding 
a fermion inverts this current, indicating that the ground state 
in this sector is at $2 ~ 27r((/3o/27r). If however $1 ~ 0, there 
is no vortivity in the Ci cycle, and the Josephson effect is the 
usual one, periodic in the flux quantum ipo = h/2e. 



{ps/2) ^{V(j)Yd?x is simply proportional to the super- 
fluid density Eh — Aps, where A is in general an 0(1) 
constant. If we have Eh < A the superconducting 
gap, then the Hopflon is the lowest energy fermionic ex- 
citation. Tunneling a single electron onto the surface 
should then spontaneously generate these phase wind- 
ings in equilibrium. Measuring the corresponding cur- 
rents (eg. via an RF squid) can be used to establish the 
presence of the Hopf texture. A daunting aspect of this 
scheme is to obtain a fully gapped superconductor with 
Ps < A- 

A different approach relies on the Josephson effect as 
illustrated in Figure [H A hollow cylinder of topologi- 
cal insulators is partially coated with superconductor on 
the top and bottom surfaces, forming a pair of Joseph- 
son junctions. A unit vortex along the Ci cycle can 
be induced by enforcing a phase difference of tt between 
the top and bottom surfaces, using the flux <I>i — ipo/2 
(where tpo — h/2e is the superconductor flux quantum). 
Now, the vorticity enclosed by the annulus determines 
the Hopf number, and hence the ground state fermion 
parity. This vorticity can be traded for magnetic flux $2 
(parameterized via /2 = 27r($2/</'o)) threading the cylin- 
der, since only Vcj) — eA is gauge invariant. Consider 
beginning in the ground state with /2 = and then tun- 
ing to f2 = 2tt. One is now in an excited state since the 
ground state at this point has odd fermion parity. This 
must be reflected in the Josephson current /. We argue 
this implies doubling the flux period of the Josephson cur- 
rent. Since / = dE{^2)/d^2, the area under the (/, $2) 
curve : d$2 /[$2] = E[ipq] - E[0] > is the excitation 
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energy which does not vanish. Hence, the Josephson cur- 
rent is not periodic in flux ipo, as in usual Josephson junc- 
tions. If we started with an odd fermion number to be- 
gin with, then the state of affairs would be reversed - the 
ground state would be achieved at multiples of /2 = 27r. 
The ground state with a particular fermion parity can 
be located by studying the slope of the current vs phase 
curve. Since the energy of the ground state increases on 
making a phase twist dl/d^2 = d^E/d^2 > 0, it is asso- 
ciated with a positive /vs<i>2 slope. This positive slope 
will be at even (odd) multiples of /2 = 2n for even (odd) 
fermion parity. If on the other hand unit vorticity was 
not induced in the cycle Ci, (eg. if $i ~ 0), the Joseph- 
son relation would be the usual one - i.e. one that is 
periodic in f2 = 2tt. This is summarized in Figure ID 

Note, a similar anomalous 47r Josephson periodicity 
was pointed out in the context pf t he 2D QSH case with 
proximate superconductivity in [2J| . We interpret this re- 
sult in terms of the fermionic nature of the solitons there 
- which lends a unified perspective. In the 3D case, the 
Hopf texture allows one to tune between the normal and 
anomalous Josephson effect by tuning the Ci vorticity 
via $1 . 



CONCLUSIONS 

The low energy field theory of the superconductor-TI 
system was derived and shown to possess a topological 
Berry phase term, which leads to fermionic Hopf soli- 
tons. We note that topological terms are particularly 
important in the presence of strong quantum fluctua- 
tions. For example, in one dimension where fluctuations 
dominate, the Berry phase term of the spin 1/2 Heisen- 
berg chainJ25|] leads to an algebraic phase. By analogy, 
it would be very interesting to study the destruction of 
superconductivity on a TI surface driven by quantum 
fluctuations. The Berry phase, or relatedly, the fact that 
a conventional insulator must break time reversal on the 
TI surface, will provide an interesting twist to the well 
know superconductor-insulator transition studied on con- 
ventional substrates 
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APPENDIX 

A. Connections to 3D non-Abelian statistics 
It is well known that vortices piercing a superconductor 
on the topological insulator surface carry Majorana zero 
modes in their cores 0, [3l ■ In the n vector representa- 
tion of Eq. dl]), this corresponds to a hedgehog defect fl5|. 
a singular configuration where the vector points radially 
outwards from the center. Although in this paper we 




(IV) 



FIG. 5: (color online) Creating a Hopf soliton by a 
hedgehog-antihedgehog pair. The blue(red) dot is the hedge- 
hog(antihedgehog). The solid and dotted lines are preimages 
of two different points on the 2-sphere. (I) creating a pair of 
hedge and antihedgehog - these could be a vortex- ant i vortex 
pair on the surface of a topological insulator (II) rotating 
the hedgehog by 2n while leaving the antihedgehog invariant 
(III) (IV) annihilating the hedgehog-antihedgehog pair. The 
final state in (IV) clearly shows linking number 1 of the two 
preimage loops, which indicates the non-trivial Hopf index. 



only work with smooth textures, we discuss below an 
indirect connection with those works. Note, we can go 
from a trivial texture to the Hopf texture by creating a 
hedgehog-antihedgehog pair, rotating one of them by an 
angle of 27r, and annihilating them to recover a smooth 
texture. This is just the Hopf texture, as can be seen in 
the Figure [5l However, as pointed out in Ref. in 
the process of rotation, the Majorana mode changes sign. 
This signals a change in fermion parity, consistent with 
our results. 

B. Calculation of Berry 's phase. We begin with: 

H = *^t>F(-iaiai)* + Af*^(ni/3o + n2P5Vx + na/Jsr?^)* 

(5) 

where, '4' is the fermion in majorana basis (8- 
component), ni is the insulator mass and n2,n3 are the 
superconductivity masses, and the standard five anti- 
commuting 4 by 4 Dirac matrices in the Majorana basis 
(where the as are symmetric and /3s are antisymmetric 
matrices): 



ai = CTz, a2 



TxCTx 



as 



-Tz(Jx 



/3o 



(6) 



We further assume the order parameters are restricted 



to unit 2-sphere: ^ ■ 



1 so that n — (711,712,713) is a 



unit vector living on S^. 

We need to show that starting from this fermionic 
model Eq.([5|) and integrating out the fermions, the ob- 
tained S'^-NLSM has an imaginary term (topological 
Berry's phase) i0H^^^g2-^(n{x,y,z,t)) with 9 = ir in the 
action (from now on we use Ht^^^m) to denote the ho- 
motopy index of a mapping.). Because this term is 
non-perturbative, in order to compute it, we need to 
embed the manifold 5^ into a larger manifold M with 
774 (M) = 0, which allows us to smoothly deform a 
H-jr^(s^){ri{x,y, z,t)) = 1 mapping to a constant map- 
ping. This means that a HTr^(^g2-^{h{x,y,z,t)) — 1 map- 
ping can be smoothly extended over the 5-dimension 
disk: V{x,y,z,t,p) : M {p ^ [0,1]) such that 

on the boundary: V{x^y,z,t,p = 1) = n{x,y,z,t) and 



7 



V{x,y, z,t, p = 0) = T/q is constant. With an extension 
V, we can perturbatively keep track of the total change 
of Berry's phase when going from a constant mapping to 
a non-trivial mapping. 

How to find a suitable M? We note the global symme- 
try of model Eq.(IS|) in the massless limit is U{l)chirai x 
SU{2)isospin, whose generators are: 

U{l)chiral ■ 75 = -ifiok SU{2)isospin ■ Vy , IbIVx , Vz} ■ 

(7) 

In our convention, /3o 7/35:75 are all anti-symmetric ma- 
trices. Starting from a given mass, for instance, (3o^ , 
one can generate the full order parameter manifold by ac- 
tionoi SUi2)isospin: ^^W^qW^^. SU[2)isospin is bro- 
ken down to U{1), the invariant group generated by rjy. 
Thus the order parameter space is SU{2)/U{1) = S^. 

Now we generalize the 8-component majorana fermion 
^ to 12-component ^. The 2-dimensional ?72:,y,z-space 
is enlarged to a 3-dimensional space, and we let the 
eight Ai-matrices (i = 1,2. ..8) of the standard SU{3) 
generators(see, e.g., page 61 in (i^]) act within this 
3-dimensional space, among which A2,A5,A7 are anti- 
symmetric while others are symmetric. And Ai_2,3 are 
just the old rjx^y^z matrices. The symmetry of the gen- 
eralized massless theory of ^': H — {—idiai)'^ is 
U{l)chirai X SU{'i)isospin, wherc the generators of the 
SU {3)isospin. are A2, A5, A7,75{Ai, A3, A4, Ag, Ag}. 

Starting from a given mass ^^/3o'4', we use 
SU{5)isospin to generate the order parameter manifold: 
*t[//3o[/t^ = ^^v^, U e SUi3),sosp^n■ It is clear 
that the 50(3) subgroup generated by A2,A5,A7 is the 
invariant group and the order parameter manifold is 
M3 = 5[/(3)/50(3). We thus embed the original or- 
der parameter manifold into M3 , and it is known that 

^4(M3) =0(0). 

The idea is to smoothly extend a H^^(^s2^^ = 1 map- 
ping over (denote by V{x,y, z,t, p)) by embedding 
5^ into Afa. In fact if we can extend a H^^(^gu(2)) = 
1 mapping over (denote by U{x,y^ z,t, p)) by em- 
bedding SU{2)isospin into SU{S)isospin, it will generate 
the extension V by V — UfioW. This is because a 
Htt4^{S^)) = 1 mapping can be thought as a combination 
of a i?^4(st/(2)) = 1 mapping and a Hopf -ff^rais^) = 1 
mapping. Such an extension U : — > SU{3)isospin has 
already been explicitly given by Wittcn (sec Eq.(9-13) in 



Ref (29! ) ■ which we will term it Witten's map. Witten's 
map was introduced to compute a i^-ff7r4(sc/(2)) topologi- 
cal Berry's phase. Basically, on the boundary dD^ — 5"* 
(p = 1), Witten's map is a Ht^^(^su(2)) = 1 mapping de- 
fined as a rotating ^^,^3(5(7(2)) — 1 soliton (by 2tt) along 
the time direction. This i?,r4(s;7(2)) — 1 mapping at p = 1 
is smoothly deformed into a trivial mapping at p = by 
embedding SU{2) into 5C/(3). 

We use Witten's map to generate V , with which we 
compute the Berry's phase perturbatively by a large mass 
expansion[l3, [l^ of Lagrangian: 

L = *t[a^ + id.,a,]^ + A/*V* = (8) 

where D = [dr + idiUi + MV] and partition function 
is Z = / V^'^V^VVer i'^^''^. After integrating out 
fermion, we obtain a NLSM of V^: L = -iTrln[ar + 
idiUi + MV]. Here the factor 1/2 is because we are in- 
tegrating out majorana fields. If there is a variation 5V , 
the variation of the imaginary part V = Im(L) is: 

Sr = ^ J dx'^e°'f^'"'TT{-f5Vdc,VdpVdjVd^VSV} (9) 

where K = f = Denoting dp ^ 84, 

after some algebra, the Berry's phase can be written in 
the fully antisymmetric way: 

r = ^ J dx^e''^'"'-'Tv{^5VdaVd^Vd^Vd^Vd,V} 

(10) 

In fact r is not fully well defined because the ambiguity 
of the extension of fi{x, y, z, t) to the 5-disk: two different 
extensions V{x,y, z,t, p) can differ by a mapping 5^ 
M3. We will soon show that this ambiguity means that 
r is well defined up to mod 2'k. 

As V is generated by [/, plugging V — U(3oU\U G 
SU{2>)isospin in Eg. (fTO|) . one can further simplify it 
by trace out the /3o,5 space. Firstly note that 
d^y = U[Wdf,U,(3o\U\ and because Wd^U is an 
element of the SU{'i)isospin Lie algebra spanned by 
A2,A5, A7,75{Ai,A3, A4,A6,A8}, [Wdf_,U,l3t)] simply picks 
out the latter five generators. After some algebra, one 
finds 
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dx^e^f"'^^''Tr{g-^d^g)^{g-^dpg)^{g-^d^g)^{g-^d^g)^{g-'d,g)^}, 



(11) 



where g is defined as the corresponding 3 by 3 5t/(3) matrix oft/: if J7 is the exponential of aiA2-l-a2A5-l-a3A7-|- 
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0475 Ai + 0575 A3 + a675 A4 + 0775 Ae + a875 As , then g is the 
exponential of ai A2 + a2X5 + a^Xr + a^Xi + a^X^ + a^Xi + 
u^Xq + agAg. g is nothing but the Witten's map. And 
{g~^df^g)± denotes the symmetric part: {g'^d^g)±^ — 
[{g-'d,g) + {g-^d,gr]/2. 

We simply need to compute F by integration. Because 
it is clear that F can only be or tt (mod 2-k), a numerical 
integration is enough to determine it unambiguously. We 
performed integration Eq. ljlip with Witten's map g by 
standard Monte Carlo approach, and find F — (1.000 ± 
0.005)7r. This proves that Hopf-skyrmion is a fermion. 
In addition it confirms that F is well-defined only up to 
mod 27r: different extensions of g can differ by a doubled 
Witten's map is known to have i?^5(sy(3)) = 1, and the 
above calculation indicate that this ambiguity only add 
an integer times 27r in F. 
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